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Abstract 

We apply the Dirac factorization method to the nonrelativistic harmonic oscillator and, more 
in general, to Hamiltonians with a generic potential. It is shown that this procedure naturally 
leads to a supersymmetric formulation of the problems under study. It is also speculated on the 
physical meaning underlying this method and it is suggested that the vacuum field fluctuations 
can be viewed as the spontaneous emission of the associated two-level system, whose quantization 
is due to the noncommuting nature of the harmonic oscillator canonical variables. 
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I. INTRODUCTION 

As stressed in Ref. [1], spin systems (i.e., any system with only two energy levels) and 
larmonic oscillators comprise two archetypes in quantum mechanics. Recent experiments 
2] are going deeper in their peculiar nature and in this paper we show that two-level spin 
systems and oscillator states are naturally entangled in a supersymmetric framework. It will 
be shown, starting from fairly simple mathematical considerations, that such a realization 
involves nothing but that the Dirac factorization method. 

The Hamiltonian of a harmonic oscillator can be written in terms of creation-annihilation 
operators according to the identity 1 

H= 1 -(p 2 + q 2 )=a+a-+ 1 - (1.1) 



where 



and 



The Pauli matrices 



in 



a ± = -L(q T tp) (1.2) 

[q,p]=i [a-,a + ] = l. (1.3) 





ox = | ; : I a 2 = | " | a 3 = ( l _° | (1.4) 

are the realization of the generators of a Clifford algebra and satisfy the identities 

[cTj, a k ] = 2 i ejkm o m , {<7j, a k } = 2 5 jk , (1.5) 

which allow us to rewrite the sum of the squares of two operators as follows 

A 2 + B 2 = (A aj + B a k ) 2 - % e jkl [A, B] a x (j ± k) . (1.6) 

This identity will be referred as Dirac factorization method. It was the breakthrough paving 



the way to the relativistic Dirac equation [3|. More recently 4j, an analogous procedure 
has been applied to get a unified view of the theory of relativistic wave equations, using 
identity ( 11. 61) as a tool to provide a different formulation of the definition of fractional oper- 
ators and derivatives. We will exploit and further expand the formalism developed in Ref. 



1 In the following we use the natural units h = c = 1 and assume, for semplicity, m — 1 and ui = 1. 



4j by applying the Dirac factorization method to obtain alternative forms of the harmonic 
oscillator Hamiltonian. This technique leads to the introduction of a set of operators which 
are not the ordinary creation-annihilation pairs and are recognized as supercharges. The 
proposed procedure brings naturally to a supersymmetric formulation of the harmonic os- 
cillator Hamiltonian and it is easily generalized to more complicated Hamiltonian forms. 
The method suggests that oscillator-like Hamiltonians are equipped, through Dirac factor- 
ization method, with a two-level structure providing a fairly transparent understanding of 
the physical role played by the algebraic grading. 

II. DIRAC FACTORIZATION AND THE HARMONIC OSCILLATOR 

To proceed in the Dirac factorization of the Hamiltonian fll.ip . let us introduce the 
following combination of matrices and differential operators 

1 ( a~\ 

E = — {qa 1 -pa 2 )=\ (2.1) 

V2 I a + / 

that, according to eqs. ( II .5p . and as a consequence of the non-commuting nature of the 
operators q and p, yields 

^ 2 = \ {(<7 2 + p 2 )-*cr l0 r 2 } , (2.2) 

and, thus, the Hamiltonian (II. ip can be rewritten as 

1 n \ 



H_ 



H = ^-^ 3 =\ 2_ 



(2.3) 



2/ 



with 

H+ = a~a + , H- = a + a~ . (2.4) 

This form of the harmonic oscillator Hamiltonian can be interpreted in terms of supersym- 



metric Quantum Mechanics (SUSY QM) 



5J because, as shown by Eq. (12.71) below, the 



operator £ can be understood as the sum of the supercharge operators associated with this 
specific problem. Let us remark that from the point of view of recovering the harmonic oscil- 
lator Hamiltonian, the expression (12. ip it's not the only one that can be adopted. However, 
it's easy to show that the general form 

S= J—^ aOj + ftgfc) (j^k), (2.5) 

y/a 2 + p 2 



can always reduced to the form (12.11) by means of a unitary transformation 

U = exp < i ^ a m o- m \ ■ (2.6) 



The term proportional to 03 in Eq. (12. 3 \ is linked to the vacuum field fluctuations. In 
the usual treatment of SUSY QM the ground state is associated only with the operator if_, 
while the present formalism suggests that the vacuum field fluctuations can be interpreted 
as a contribution emerging from a kind of population inversion. We note indeed that, apart 
from its mathematical role, the operator S is amenable for a transparent interpretation in 
physical terms. The joint use of the properties of Pauli matrices and ladder operators allows 
us to cast Eq. ( 12.1 j) in the form 

S = a~ a + + a + a_ (2.7) 

with 

°+ = o( flr i + io "2) =( , o- = 0(0-1 -icr 2 ) = {2A 





The operator ( 1 2 .7p is formally equivalent to the interaction potential in Jaynes-Cummings 
Hamiltonian [6], describing the interaction of a quantized field with a two- level system, 
whose energies differ by the characteristic gap of the harmonic oscillator spectrum 2 . We can 
therefore say that the use of the Dirac factorization method and the quantized nature of 
the Hamiltonian operator induces, in quite a natural way, an analogous quantized two level- 
structure. This, in turn, can be interpreted as the physical origin of the graded algebraic 
nature of the model we are developing. 

The algebraic structure, underlying the Hamiltonian (I2.3p . requires the embedding of 
bosonic (a*) and fermionic (Pauli matrices) operators which close up to form a super- 
algebra. Following Refs. [7H9J], we introduce the operators 

(2.9) 



U± = a a± 


v± 


= a T o± 




that satisfy the commutation relations 








[U + ,U.]=H<T 3 -l 




\y+,v-\ 


s 2 


[U±,V±]=T°l 




[U±,Vr] 


= ±2K ± a 3 



(2.10) 



This analogy is an a posteriori justification for our particular choice (|2.f [) for the operator E. 



where we have introduced the two operators 

K ± = 1 -(a ± )\ (2.11) 

that, together the Hamiltonian operator, generate the SU(1,1) algebra: 

[H,K ± ] = ±2K ± [K + ,K_] = -H (2.12) 

and allows us to recover the full structure of the ortho-symplectic algebra osp(l|2) with the 
remaining commutation brackets (the not mentioned brackets are zero) 

[ K± U t } = tV t [K t ,V ± } = tU ± . (2.13) 

As anticipated, the operators V± play the role of supercharges of the system. Indeed, 
from Eq. ( 12.71) . we can write £ = V+ + V-, and, since V± = 0, the supersymmetric part of 
the Hamiltonian is given by S 2 = {V+, V-}. 

By using the previous identities, it's easy to show that the Heisenberg equation of motion 
for the operators a^ are: 

a 1 = & 2 = & 3 = 2{q,p} = 4i(K + -K_) (2.14) 

from which it's easy to show that 

03(f) oc sin 2 (V2t) , (2.15) 

that represents the vacuum field spontaneous emission. As for the time evolution of the 
operator £ we obtain 

£ = -£ (2.16) 

i.e., a simple harmonic oscillation, physically associated with the emission and the absorption 
of a photon. 

We have so far shown that the Dirac factorization method dresses the harmonic oscillator 
with a spin- like (two- level structure), yielding a graded Lie algebraic structure. It is therefore 
natural to ask about the real physical meaning of such a superimposed structure. The 
main representative of the fermionic structure is the vector a = (ax, o-i, 03), which may be 
considered as rather artificial. However, apart from the fact that it emerges in a very natural 
way from the mathematical procedure, its role is physically understandable and should be 



interpreted as that of a quantized two-level structure induced by the noncommuting nature 
of the q and p variables. More in general, we can view the vacuum as an ensemble of coupled 
two-level systems continuously emitting and absorbing a photon. 

Before closing the section, we briefly show how the Dirac factorization proceeds in the 
case of the slightly more complicated Hamiltonian 



H=~ [ P 2 + f(q)} ■ 
If, in analogy with Eq. ( 12. ip . we introduce the operator 



(2.17) 



T 



V2 



(vTS) 



«7i ~pcr 2/ 




with 



we can write 



H = T 2 



Also in this case 



1 /'(g) 



a ± = -^(v75)t<p) 



(A- A-- 1 - ™ 



(2.18) 



(2.19) 



o% 



\ 



4 v / 7¥ 

A + A- 



\ 

1 f'(q) 

4 v75J/ 



(2.20) 



V_ 



(2.21) 



A- \ / 

0/ \ A + 

are recognized as super-charges associated with the Hamiltonian ( j2.17p . and the quantities 



f± = f(q) ± 7 



(2.22) 



4 v / 7M 

are interpreted as super-partners potentials. The example which follows has a twofold mo- 
tivation, i.e., an explicit application of the techniques we have developed and the use of 
nonstandard special functions emerging from the analysis of the problem under study. We 
consider the case of the quartic oscillator, i.e. f(q) = Ag 4 , for which the ladder operators 
are 

A ± = ±=(V\q 2 T ip) (2.23) 

with 

[A~,A + ]=2V\q. (2.24) 
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The states can be denned as follows 

^ = -L { A + ) n |0> (2.25) 

Vn! 

where, from the condition A~ |0) = 0, for the vacuum it is possibile to deduce the following 
expression 

|0)cxexp(-^g 3 j (q>0). (2.26) 

The explicit form of the functions f)2.25|) can be obtained by using the identities (see Ap- 
pendix) 

( A+ r = ^E (t) (-i)**& (^« 2 . ^«. x) 9 « (2 - 27) 

a j e -/Ir'/3 = ff (3) | _VA 9 2 , -VX«, -^ J c-^«" /3 (2.28) 

where i?A are the third-order Hermite polynomials whose definition is given in Appendix. 
The use of the addition theorem 



5Z ( jfc ) H n-k( x u x ^ x s) HJ®(yi, y 2 , y 3 ) = H®(x t + y 1} x 2 + y 2 , x 3 + y 3 ) 

fc=0 ^ ' 

and taking into account that (— l) fc H\, \x, y, z) = H\, \—x, y, —z), finally yields 



(2.29) 



VnOC-^H^flVX^O^Vxj e-^ 3 /3. (2 . 30) 

The functions (p n are significantly different from the ordinary harmonic oscillator functions, 
and the discussion of their properties is out of the scope of this paper. 

In the forthcoming section we will consider more physical examples yielding further ele- 
ments supporting the physical reality of the mathematical devices discussed so far. 

III. PHYSICAL EXAMPLES 

We have mentioned that the quantized nature of the spin-like structure introduced in the 
previous sections reflects the quantum nature of q and p variables. The same happens with 
the so-called Landau states, emerging in the quantum analysis of the motion of a particle 
with electric charge e in a classical magnetic field of intensity B |l| . 



By identifying the z-axis of the reference frame with the direction of the field, the asso- 
ciated vector potential writes A = (0,Bx,0). This choice is not unique but it ensures that 
the magnetic field is orthogonal to the plane of motion, i.e., p = (p x ,P y , 0). The relativistic 
Hamiltonian operator ruling such a process can be written as 3 



H - p-eX) +(mc) 2 



= c^Jp 2 x +p 2 y -2eBxp y + (eBx) 2 + (mc) 2 . (3.1) 

Since [p y ,H] = 0, the operator p y can be replaced by its eigenvalue h k y , and the Hamiltonian 
can be rewritten as 



H = c Jp 2 + m 2 u 2 X 2 + {h ky) 2 — mu 2 x 2 B + (m c) 2 . (3.2) 

with 

lei B hk v 

uj c = , x B = - , X = x - x B ■ 

mc mui c 

By introducing the operator 



W = p x (J\ + muj c X a 2 + \/(hk y ) 2 — mcu 2 x 2 B + (mc) 2 a s (3.3) 

it's easy to show that 

H = c^W 2 -mhu c a- i . (3.4) 

In the common experimental situations, the term proportional to 03 can be neglected, since, 
for values of the magnetic field around 0.1 T (or even larger) it is of the order of meV, while 
the term c W varies in the region of hundreds of keV. Under this approximation, the study 
of the relativistic Landau states reduces to a Jaynes-Cummings problem. The Hamiltonian 
(13. 4p describes the dynamics of a two-level system with a level spacing fixed by the strength 
of the magnetic field. The supersymmetric nature of the Jaynes-Cummings model has been 
discussed elsewhere [6( and we will not dwell on it. Here we note that a physical realization 



of such a system is given, for example, by a Free Electron Laser (F 
relativistic beam of electrons propagates inside an axial magnetic field 



L) source where a 



IOJ. The possibility 



of using such a point of view to construct a laser-like theory for the FEL-like devices has 



been partially considered in [ll| , and will be the topic of a future speculation. 



In this section we restore the mass m of the particle and the constants H and c. 



!4 



Before closing this section we will reconsider a point partially touched in Ref. |4|, and 
in previous sections, concerning the evolution of the spin-like system associated with the 
harmonic oscillator via the Dirac factorization method. To simplify the problem we will 
refer to the relativistic 1-dimensional Hamiltonian 

H = c y /p 2 + (mc) 2 (3.5) 

that, for example, can be factor ized as follows 

H = c{pGi + mc<7 3 ) (3.6) 



i, 



The physical meaning of the above Hamiltonian has been discussed in Ref. j4|, where it has 
been stressed that it should not be confused with the Dirac Hamiltonian nor with the Pauli 
counterpart. Here, we will use it as a toy model to get a further support to the previous 
speculations. 

The equations of motion for the vector a are 

— a = tixa n = -Qoc,0,mc 2 ). (3.7) 

This equation describes a purely quantum motion (with no classical counterpart) which 



should be understood as a kind of zitterbewegung |4|, [12] (a trembling motion due to the 
interference between negative and positive states contained in the Hamiltonian (13. 6p ). Even 
though there is not any explicit presence of ladder operators, the Hamiltonian (I3.6P contains 
a hidden two-level structure and a supersymmetric underlying algebra. In fact, by expressing 
the momentum in terms of the ladder operators, we get (g = c/y2h) 

H = mc 2 cr 3 + ihg (a + + a") (<r + — er_) , (3.8) 

which writes as a Jaynes-Cummings Hamiltonian, but without the rotating wave approxi- 
mation assumption. 

IV. CONCLUDING REMARKS 

As it is shown in sec. [TTJ the method of factorization exhibits, in a fairly natural way, 
all the essential features of SUSY QM: i) the vacuum field energy is factorized out; ii) two 
potentials with isospectral properties are recovered. But, as often happens, nothing is really 

9 



new and the procedures leading to super-potentials traces back to methods known well before 
to the birth of quantum mechanics itself |13l . |l4j 4 . In order to appreciate this point, let us 
consider the following second-order differential equation: 

z"(x)+fi_(x)z(x) =0. (4.1) 

If the "potential" /i_(x) is expressible as 

^{x) = -\cf z {x)- l -${x) ) (4.2) 

i.e., is solution of a Riccati equation, the solution of Eq. ( 14. ip can written as 

z(x) =exp j^ /V<K0} u-(s) (4.3) 

with the function u_(x) satisfying the differential equation 

u'L(x) + (j)(x) u'_(x) = 0. (4.4) 

Therefore we get 

z(x)=expl- I d^0(O} / dr/expj-/ ducf>{u)\ . (4.5) 

Furthermore once 4>(x) is fixed, we can define a second "potential" 

^ + (x) = -^ 2 (x) + ^'W (4-6) 

(the super-partner of /i_(x), according to the present terminology) which specifies the solu- 
tions of a second differential equation, obtained from the first by just replacing 0(x) with 
-<f>{x). 

The outlined procedure is essentially the Liouville method to reduce a second-order dif- 
ferential equation to its standard form 5 . However, it is remarkable that all the features 



4 In particular, in Ref. [13| the problem has been treated by showing that a second-order differential 
equation with non-constant coefficients can be written in terms of bi-orthogonal partner solutions, a 
posteriori recognized as supersymmetric components. 

5 We remind that, given the differential equation y" +a(x) y'(x) + b(x) y(x) = 0, the Liouville transformation 



y(x)=cxpl-- d£a(m v(a 



r 



educes the initial equation to v"{x) + c(x) v(x) = 0, where c(x) = b(x) — j [a 2 (x) + 2 a'(x)]. 
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concerning the supersymmetry are recovered in a very natural way by just applying the 
Dirac factorization method, which leads to the two-level Hamiltonian (I2.17P with the role of 
the vacuum field factorized out from the very beginning. 

Let us come back to the operator T defined in Eq. (12.1 81) . Also in this case, we can write 

T = A- o+ + A + a_ , (4.7) 

and, therefore, it can be be interpreted as the interaction between a two-level system and 
the bosonic field defined by ladder operators given in Eq. (I2.19p . 

More in general, we can conjecture that the Jaynes-Cummings model can be generalized 
by the following interaction Hamiltonian 

# JC = T - 1 U {q) a, L(q) ex -^jL j . (4.8) 

As a consequence of the more complicated commutation relations involved, the algebraic 
nature of this operator is less direct than in the case of the ordinary J-C operator. However, 
the study of quantum states ruled by the interaction Hamiltonian (14. 8 p can be done using 
the evolution operator associated with the Hamiltonian The properties of this operator can 
be discussed by using the methods developed in the past to treat evolution problems in 



quantum mechanics (see Refs. [15|, |16(). For example, using the so-called symmetric-split 
decomposition method [17] . the evolution operator can be approximated as follows 

N 

*(t)=n^*(°) ( 4 -°) 

with 

Uj = e~^ T/2 e-^ (9)<r3 e~ itj T/2 + 0{t)) (4.10) 

The use of standard identities for exponential operator (see Ref. [161]). along with the iterated 
application of the evolution operator in Eq. (14. 9p . yields an efficient way of calculating the 
evolution of these quantum states. 

In this paper we have touched different topics, whose underlying leitmotif is the Dirac 
factorization. We have discussed the relativistic quantum mechanics, the zitterbewegung, 
the Jaynes-Cummings model, the relativistic Landau levels, and the Dirac oscillator (even 
though not explicitly mentioned). Most of these effects are difficult to be subjected to exper- 
imental investigation. For example, the electron zitterbewegung exhibits oscillations at very 

11 



large frequencies (~ 10 21 Hz), not accessible to currently available experimental techniques. 



Such drawback is overcome in the context referred to as quantum simulation 18 



where the 



trembling motion may occur, for example, in crystalline solids when their band structure is 



represented by a two-band model reminiscent of the one- dimensional Dirac equation 19 ]. 



The explicit realization of the simulation 18|, |20| has provided strong indications for the ex 



istence of such a genuine quantum behavior. Further experiments using quantum simulator 



techniques are suggested to test other quantum paradoxes like the Klein paradox 21]. 

We believe that the topics treated in this paper can be framed within the context of the 
quantum phenomenology which can be experimentally tested using quantum simulators. A 
more specific analysis in this direction will be developed elsewhere. 

Appendix A 

In section IIVI we have used the higher-order Hermite polynomials to express the ordered 
form of the operator O n = (d x + P(x)) n . These polynomials are defined through the gener- 
ating function 

^-^(x 1 ,x 2 ,---,x m )=exp^x fc tH . (Al) 

n=0 ' I k=l J 

and can be constructed recursively according to formula 

[n/m] k 

X, 



H<THx 1 ,x 2 ,~-,x m ) = n\ J2 in Xm iAi ^fa^-.^i)- (A2) 

fc=0 K -\ n ~ mK )- 

By using the generating function methods, we can define the following operator 

00 j-n 

E(x, t) = Y j -O n = exp{t [d x + P(x)]} (A3) 

n=0 

that satisfies the differential equation 

d t E{x, t) = [d x + P(x)} E(x, t) E{x, 0) = 1 (A4) 

whose solution is 

E(x,t) =expltd x + I dt'P{x-t')\ = exp | f dt' P(x + t - t')\ e tdx . (A5) 

In the case P(x) = a x 2 , performing the integral and using the Eq. (I A II) . we easily obtain 

OO f. 

E(x,t) = J2^H^(ax\ax^) e td * (A6) 



12 



and, thus, 

a.=YV n W>. (a^.ax. 






H® k (ax*,ax^)d x k . (A7) 



fc=0 

As for the successive derivatives of the function e~ x , by applying the generating function 
method we get 



00 + n °° j-n 

t- ^^, ~3 /a _™3 _f™_l_+'\3 X *> l> 



Y'-d" e~ x3 = e td * e~ x3 = e"^ 3 = V - H®(-3 x\ -3 x, -1) e"* 8 (A8) 

^— ' n\ ^— ' n! 

n=0 n 

and, therefore, from the comparison of the same powers of t, the second identity in (12.281) . 
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